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We study the thermodynamic stability at low temperatures of a series of alkali metal/zinc double- 
cation borohydrides, including LiZn(BH4)3, LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5, KZn(BH4)3, 
and KZn2(BH4)5. While LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5 and KZn(BH4)3 were recently 
synthesized, LiZn(BH4)3 and KZn2(BH4)5 are hypothetical compounds. Using the minima-hopping 
method, we discover two new lowest-energy structures for NaZn(BH4)3 and KZn2(BH4)5 which 
belong to the C2/c and P2 space groups, respectively. These structures are predicted to be both 
thermodynamically stable and dynamically stable, implying that their existence may be possible. 
On the other hand, the lowest-energy PI structure of LiZn(BH4)3 is predicted to be unstable, 
suggesting a possible reason elucidating why this compound has not been experimentally identified. 
In exploring the low-energy structures of these compounds, we find that their energetic ordering is 
sensitive to the inclusion of the van der Waals interactions. We also find that a proper treatment 
of these interactions, e.g., as given by a non-local density functional such as vdW-DF2, is necessary 
to address the stability of the low-energy structures of these compounds. 

PACS numbers: 61.66.-f, 63.20.dk, 61.05.cp 
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I. INTRODUCTION 

Hydrogen is a promising alternative to fossil fuels as 
it can provide a high density of clean energy. To use 
hydrogen as a fuel, a safe and efficient storage method 
is crucial, especially for mobile applications.^^ Many 
solid state materials, including complex borohydrides 
M(BH4)^, have been studied for hydrogen- storage pur- 
poses (here M is a metal of valence m). However, none of 
these compounds currently meets all the requirements for 
a hydrogen storage medium. For examples, the dehydro- 
genation of alkaline-earth metal borohydrides M(BH4)2 
is generally too slow, while alkali metal borohydrides 
MBH4 are thermodynamically too stable. Attempts to 
improve their performance, e.g., by using additives, by 
partial substitutions, or by subjecting them to confine- 
ment at the nanoscale, have been met by limited success.^ 

A large number of double-cation borohydrides 
MiNj{B'¥{^)mi+nj were recently synthesized and pro- 
posed as candidates of hydrogen-storage materials 
[N is also a metal of valence n). Several exam- 
ples of such compounds include Li/K borohydrides,^ 
Li/Ca borohydrides,^ Li/Sc borohydrides,^ Na/Sc 
borohydrides,^ K/Sc borohydrides,^^ Na/Al 

borohydrides,-^-^ K/Mn and K/Mg borohydrides,^^ K/Y 
borohydrides,^^ Li/Zn and Na/Zn borohydrides, ^^^^ 
K/Zn borohydrides.^^ These compounds, whose hydro- 
gen storage capacity is relatively high, were suggested 
to have decomposition temperatures lying somewhere 
between those of the corresponding single-cation 
borohydrides. ^'^'^'^^'^^'^^'^^ Such compelling feature 
would allow for the adjustment of the decomposition 



temperature by selecting an appropriate cation combi- 
nation. Theoretical studies on these newly developed 
materials were subsequently carried out.^^'^^^^ 

In this work, we consider a series of compounds of 
Li/Zn, Na/Zn, and K/Zn double-cation borohydrides 
of which LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5 and 
KZn(BH4)3 were experimentally synthesized. ^^^^ At 
ambient conditions, LiZn2(BH4)5 crystallizes in an or- 
thorhombic Cmca (no. 64) phase as first determined by 
powder x-ray diffraction (XRD) analysis^^ and then re- 
fined by powder neutron diffraction, Raman and NMR 
spectroscopy.^^'^^ Next, NaZn(BH4)3 and NaZn2(BH4)5 
were found to be in two phases which both belong to the 
monoclinic P2i/c space group (no. 14).^^'^^ At lOOK, 
KZn(BH4)3 was determined to crystallize in a rhombo- 
hedral R3 phase (no. 146).^^ Two related hypotheti- 
cal compounds that are also considered in this work are 
LiZn(BH4)3 and KZn2(BH4)5. While both of them have 
not experimentally been observed yet, a P2i/c phase (no. 
14) and a PI (no. 1) phase of LiZn(BH4)3 were theo- 
retically proposed in Refs. 15 and 18, respectively. No 
information on KZn2(BH4)5 is available in the literature. 

Some of the structural phases experimentally pro- 
posed for these compounds were recently suggested to 
be thermodynamically unstable at OK by a theoretical 
study based on ab initio calculations.^^ In particular, 
the unrefined Cmca phase of LiZn2(BH4)5 proposed by 
Ref. 14 was suggested to be slightly unstable with re- 
spect to the decomposition into Zn(BH4)2 and the PI 
phase of LiZn(BH4)3. Furthermore, the P2i/c phase of 
NaZn(BH4)3 was found^^ to be unstable with respect to 
the decomposition into NaBH4 + Zn(BH4)2. However, 



it turns out that the correct stabihty of LiZn2(BH4)5 
has to be examined by using the Cmca phase which 
was refined by Refs. 15,16. Also, because the refined 
Cmca phase of LiZn2(BH4)5 was noted^^ to be lower 
in energy than the unrefined Cmca phase, the conclu- 
sion of Ref. 18 that the PI structure of LiZn(BH4)3 
is stable has to be re-examined. In addition, the re- 
cent proposal of the PI phase of NaZn(BH4)3,^^ which is 
more stable than the experimentally synthesized P2i/c 
phase, ^^ suggests that unexplored low-energy structures 
of the Li/Zn and Na/Zn borohydrides may exist. Fi- 
nally, as KZn(BH4)3 was synthesized very recently,^^ it 
is worth to explore the (additional) possible low-energy 
structures of KZn(BH4)3 and its related compounds, e.g., 
KZn2(BH4)5. 

These above open issues on the thermodynamic sta- 
bility of the Li/Zn, Na/Zn, and K/Zn borohydrides are 
examined in this work. In order to discuss the sta- 
bility of these compounds, we search and discover a 
large number of low-energy structures using the minima- 
hopping method. ^^'^^ Van der Waals (vdW) interactions, 
which were recently shown^^ to be important in de- 
scribing the stabilization of the low-energy structures of 
Mg(BH4)2, is properly considered using vdW density- 
functional theory. ^^^^ Finally, we study the thermody- 
namical stability of these compounds at finite temper- 
atures by determining the vibrational free energy from 
the phonon frequency spectra of the examined low-energy 
structures within the harmonic approximation. 



II. METHODS 

First-principles calculations reported in this work 
were performed at the level of the density functional 
theory (DFT)^^'^^ within the Perdew-Burke-Ernzerhof 
(PBE) generalized gradient approximation.^^ The pro- 
jected augmented wave formalism^^ as implemented in 
Vienna Ah Initio Simulation Package (vASP)^^^^ was 
employed. The valence electron configurations of Zn, B, 
Li, Na, and K used for calculations are 3(i^^4s^, 2s^2p^ ^ 
ls^2s^, 25^2p^35^ and 35^3p^45^, respectively. The con- 
vergence of the DFT energy ^^ft ^^^ ensured by a 
Monkhorst-Pack k-point mesh^'^ from 5x5x5 to 9x9x9, 
depending on the simulation cell sizes, for sampling the 
Brillouin zone and a kinetic energy plane wave cutoff 
of 900 eV. Atomic and cell variables were simultane- 
ously relaxed until the residual forces were smaller than 
0.01 eV/A. The effect of vdW interactions of all the ex- 
amined structures was explored in separate calculations 
using the modified Langreth-Lundqvist non-local den- 
sity functional (vdW-DF2)^^^^ as implemented in the 
PWSCF code of the quantum espresso distribution.^^ 
The rPW86 pseudopotentials^^'^^ of the PsLibrary at an 
energy cutoff of 80 Ry were used in the calculations for 
the DFT energy £^dft' taking into account the vdW in- 
teractions. The space groups of the examined structures 
were determined by findsym^^ while some figures were 



prepared by VESTA. ^^ 

We tested our computational scheme by opti- 
mizing the experimentally reported structures of 
LiZn2(BH4)5,^^"^^ NaZn(BH4)3/'^'^^ NaZn2(BH4)5,^'^'^^ 
KZn(BH4)3,^^ LiBH4,'^^ NaBH4,'^'^ and KBH4.'^^''^^ The 
optimized structures are summarized in Table I, show- 
ing that the Cmica structure of LiZn2(BH4)5 obtained 
from the PBE optimization agrees very well with ex- 
perimental data.^^ Moreover, our results for the P2i/c 
phases of NaZn(BH4)3 and NaZn2(BH4)5 are consistent 
with those reported in Ref. 18, which are generally in 
good agreement with the experimental values, ^^ except 
for some discrepancies on h of NaZn(BH4)3, and on a of 
NaZn2(BH4)5. While the origin of these discrepancies is 
still unclear, ^^ we note that a similar discrepancy on h 
of the P2i/c phase of NaZn(BH4)3 was already reported 
in Ref. 18. In addition, c calculated for the R?> phase 
of KZn(BH4)3 is about 10% larger than the experimen- 
tal value. ^^ Presumably, this deviation may come from 
the finite temperature (lOOK) at which the KZn(BH4)3 
samples were analyzed and/or the fact that a, 6, and 
c were all fixed in the DFT calculations subsequently 
performed to confirm the space group i?3 determined 
experimentally. ^^'"^^ For LiBH4, NaBIl4, and KBH4, the 
calculated results agree very well with the experimen- 
tally measured values. With vdW-DF2, as also shown 
by Table I, the agreement between the calculated re- 
sults and the experimental data is improved. Similar to 
the PBE results, the lattice parameters of LiZn2(BH4)5, 
LiBH4, NaBH4, and KBH4 optimized with vdW-DF2 
agree well with the experimental values. Furthermore, 
the discrepancies previously observed for NaZn(BH4)3, 
NaZn2(BH4)5, and KZn(BH4)3 with the PBE optimiza- 
tion are strongly diminished. These observations clearly 
show that the long-ranged vdW interactions in the ex- 
amined structures can be well captured by the non-local 
vdW-DF2 functional used in this work. 

Unconstrained searches for low-energy structures of 
the examined compounds were performed by the minima- 
hopping method {MRM)?^^'^^ In this method, the PBE 
energy landscapes are explored by consecutive short 
molecular dynamics steps followed by local geometry re- 
laxations. The initial velocities for the molecular dy- 
namics runs are chosen approximately along soft mode 
directions, allowing efficient escapes from local minima, 
and aiming towards the global minimum. This method 
was successfully applied in a wide range of material struc- 
ture predictions. ^^^^ Some of the theoretical predictions 
by MHM on neutral Si clusters and four-fold defects in 
silicon were recently confirmed by experiments. ^^'^^ 

Phonon frequency spectra of the examined struc- 
tures were calculated with the PBE functional using the 
PHONOPY package. ^^'^^ Given a sufficiently large relaxed 
super cell, finite atomic displacements with an amplitude 
of 0.01 A were introduced. The atomic forces within the 
super cells were calculated with VASP and the phonon 
frequencies were then calculated from the dynamical ma- 
trix, given in terms of the force const ants. ^^ The longi- 



TABLE I: Geometrical parameters (lattice parameters a, b, and c in Angstrom, angle /3 in degree, and volume difference AV in 
percent) obtained by computationally optimizing the experimentally-reported low-energy structures of the examined compounds 
and several related compounds. References are given for the experimental data, which are presented for the comparison purpose. 



Compound 


Space group 




Calculations with PBE 


Calculations with vdW-DF2 




Experiments 




a(A) 


b{A) c(A) 


PC) AV{%) 


a(A) 6(A) c(A) PC) AV{%) 


a(A) 


6(A) c(A) PC) 


Ref. 


LiZn2(BH4)5 


Cmca 


8.63 


18.03 15.43 


90 


0.9 


8.76 17.62 15.78 90 


2.4 


8.62 


17.90 15.41 90 


[161 


NaZn(BH4)3 


P2i/c 


8.02 


5.38 18.57 


101.5 


14.3 


8.84 4.65 17.14 100.0 


1.1 


8.27 


4.52 18.76 101.7 


[141 


NaZn2(BH4)5 


P2i/c 


10.45 


16.19 9.06 


112.4 


6.4 


9.07 17.05 9.36 112.2 


1.7 


9.40 


16.64 9.14 112.7 


[141 


KZn(BH4)3 


R3 


7.80 


7.80 12.24 


90 


16.5 


7.78 7.78 11.88 90 


12.6 


7.63 


7.63 10.98 90 


[171 


LiBH4 


Pnma 


7.30 


4.39 6.60 


90 


-2.4 


6.96 4.54 6.72 90 


-2.1 


7.18 


4.44 6.80 90 


[431 


NaBH4 


P42/nmc 


4.34 


4.34 5.88 


90 


0.6 


4.36 4.36 5.86 90 


1.2 


4.33 


4.33 5.87 90 


[441 


KBH4 


P42/nmc 


4.76 


4.76 6.68 


90 


0.6 


4.72 4.72 6.60 90 


-3.4 


4.77 


4.77 6.69 90 


[451 



tudinal optical/transverse optical splitting was not con- 
sidered because its effects were reported to be negligible 
for a wide variety of hydrides. ^^'^^ 

For a double-cation borohydride, we computed the for- 
mation energy to examine if it is stable with respect to 
the decomposition into the corresponding single-cation 
borohydrides and related compounds. As an example, 
the formation energy of LiZn2(BH4)5 was determined by 



AF =F[LiZn2(BH4)5] 

-{xFlUBR^l + (1 - x)F[Zn(BH4)2]} 



(1) 



where x = 2/3 is the Zn cation composition of 
LiZn2(BH4)5. Here, F[LiZn2(BH4)5], F[LiBH4], and 
F[Zn(BH4)2] are the free energies of one formula unit of 
LiZn2(BH4)5, LiBH4, and Zn(BH4)2, respectively. The 
examined structure of LiZn2(BH4)5 is stable if AF is 
negative; otherwise, it would decompose into LiBH4 and 
Zn(BH4)2. At finite temperature T, the vibrational free 
energies were determined within the harmonic approx- 
imation from the phonon frequency spectra of the ex- 
amined structures. The formation energy was also cal- 
culated from F^DFT' i-^-5 ^^^ DFT energies of the struc- 
tures, taken into account the vdW interactions. To de- 
termine the formation energies AF according to Eq. 1, 
we used the /4i22 phase (no. 98) of Zn(BH4)2/^'^'^'^^'^^ 
the Pnma phase (no. 62) of LiBH4,^^ and the PA2/mnc 
phase (no. 137) of both NaBH4 and KBH4>'^''^^'^^ 



III. LOW-TEMPERATURE STRUCTURES OF 

ALKALI METAL/ZINC DOUBLE-CATION 

BOROHYDRIDES 

We searched for low-energy structures of LiZn(BH4)3, 

NaZn(BH4)3, KZn(BH4)3, LiZn2(BH4)5, NaZn2(BH4)5, 
and KZn2(BH4)5 by using the minima-hopping method. 
For the first three compounds, we performed searches for 
1 f.u. (17 atoms) and 2 f.u. (34 atoms) while for the 
last three compounds, searches were performed with 1 
f.u. (28 atoms). Starting from random input structures, 
we discovered a large number of new low-energy struc- 
tures for these compounds. Additional searches were 



also performed by relaxing several lowest-energy struc- 
tures discovered for each compound after replacing their 
cations by the cations of other appropriate species. Our 
MHM runs were performed with the PBE functional and 
the structures discovered were then examined using the 
vdW-DF2 functional. 

We show in Table II the densities and the formation 
energies, determined with PBE and vdW-DF2, of all 
the structures discovered together with those already re- 
ported in the literature. We note that our results for 
AF^YT ^^^ ^^^ ^^ phase of LiZn(BH4)3 and the unre- 
fined Cmca phase of LiZn2(BH4)5 are consistent with 
those reported by Ref. 18. The crystallographic infor- 
mation of all the structures we discovered can be found 
in the Supplemental Material. ^^ Simulated XRD patterns 
performed by FULLPROF^^ are also shown in the Supple- 
mental Material, ^^ indicating that the low-energy struc- 
tures reported in this work are different from those al- 
ready reported in the literature. 



A. Lithium/zinc borohydrides 

We optimized the refined Cmica structure of 
LiZn2(BH4)5 and found that this structure is thermody- 




Cmca structure 



l-4m2 structure 



FIG. 1: (Color online) Cmca (left) and /4m2 (right) phases 
of LiZn2(BH4)5. Lithium and zinc atoms are shown by green 
and gray spheres while [BH4]~ complex anions are shown 
by green tetrahedra. Zn-Zn interlinked "bonds" are used as 
guides for the eyes. 



TABLE II: Summary of the available structures of alkali metal/zinc double cation borohydrides. AF^p^ , AFq^^ , AF^^q^, 
and Fdf^? given in kJ/mol cation, are the enegies of formation determined from ^dft, ^ok, ^iook, and ^dft? respectively. 
These structures are indicated by their space groups. Densities p are given in g/cm^. References are shown for structures not 
originating from this work. 



Compounds 


Space group (no 


.) p(g/cm^) 




A T?^^^ 


A T?^^^ 
^^lOOK 


A TTlvdW 


References 


LiZn(BH4)3 


Pl(l) 


0.80 


-6.70 


-6.89 


-7.45 


-4.95 


[18] 


LiZn2(BH4)5 


Cmca (64) 


1.18 


-10.14 


-10.67 


-10.61 


-19.05 


[15,16] 




/4m2 (119) 


0.57 


-11.70 


-12.96 


-13.59 


-1.08 




NaZn(BH4)3 


P2i/c (14) 


1.28 


-6.55 


-8.43 


-8.74 


-4.90 


[14] 




Pl(l) 


0.91 


-7.83 


-9.25 


-9.92 


-4.29 


[18] 




C2/c (15) 


1.00 


-10.90 


-11.76 


-12.05 


-9.56 




NaZn2(BH4)5 


P2i/c (14) 


1.15 


-12.08 


-13.91 


-14.58 


-11.47 


[16] 




/4m2 (119) 


0.54 


-11.87 


-13.58 


-14.87 


1.90 




KZn(BH4)3 


R3 (146) 


1.15 


-7.31 


-9.28 


-10.21 


-8.86 


[17] 




P63/m (176) 


1.05 


-11.08 


-12.60 


-13.32 


-5.10 




KZn2(BH4)5 


P2(3) 


1.09 


-8.84 


-10.75 


-11.39 


-10.97 





namically more stable than the unrefined Cmca structure 
by :^ 7.8 kJ/mol cation. Additionally, from the MHM 
searches for LiZn2(BH4)5, we discovered a new tetrago- 
nal IAm2 structure (no. 119), which is more stable than 
the refined Cmca structure by c^ 1.5 kJ/mol cation. Fur- 
ther calculations with the PEE functional indicated that 
/477i2 is the thermodynamically most stable structure of 
LiZn2(BH4)5 at temperatures up to lOOK (see Table II). 
While the Cmtca structure consists of two identical inter- 
penetrated three-dimensional (3D) frameworks, ^^^^ the 
IAmi2 structure consists one of these two frameworks with 
differences in the sequence of Zn atoms and Li atoms. 
Therefore, it is not surprising that the /4m2 structures 
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FIG. 2: (Color online) Density of phonon states of (a) the 
/4m2 structure of LiZn2 (6114)5 and (b) the PI structure of 
LiZn(BH4)3. 



has a very low density with p = 0.57 g x cm~^, about 
one half of the density p = 1.18 g x cm~^ of the Cmtca 
structure. An illustration for the Cmica and IAmi2 struc- 
tures is given in Fig. 1, clearly showing the difference in 
the geometry and density of the Cm^ca and /4m2 struc- 
tures. The density of phonon states of the IAmi2 struc- 
ture is plotted in Fig. 2(a), indicating that this struc- 
ture is dynamically stable because no phonon mode with 
imaginary frequencies was observed in the whole Bril- 
louin zone. 

For the hypothetical compound LiZn(BH4)3, we opti- 
mized the P2i/c and PI structures which were theoreti- 
cally proposed in Refs. 15 and 18, respectively. We found 
that for the PI structure, AP^p|^ c^ —6.7 kJ/mol cation 
while for the P2i/c structure, AP^ff 2^ 19 kJ/mol 
cation. The positive formation energy of the P2i/c struc- 
ture hints that this structure is unstable and is an in- 
correct structural model. ^^ On the other hand, PI is the 
lowest-energy structure of LiZn(BH4)3, as also confirmed 
by the MHM runs we performed for this hypothetical 
compound. The PI structure is dynamically stable, as 
demonstrated by its density of phonon states shown in 
Fig. 2(b). 

We then examined the thermodynamic stability of the 
PI structure of LiZn(BH4)3 and the Cmtca and the IAmi2 
structures of LiZn2(BH4)5 based on the convex hull con- 
structed from the the formation energies of the thermo- 
dynamically most stable structures calculated at OK and 
lOOK with the PBE functional. Fig. 3(a) shows that in 
agreement with the conclusion of Ref. 18, the unrefined 
Cmtca phase is unstable with respect to the decomposi- 
tion into LiZn(BH4)3 and Zn(BH4)2. Because the refined 
Cmtca and IAmt2 structures of LiZn2(BH4)5 are well lower 
in energy than the unrefined Cmtca phase, they remain 
thermodynamically stable at temperatures up to lOOK 
[see Figs. 3(a)-(c)]. The PI structure of LiZn(BH4)3, 
on the other hand, is always unstable and would decom- 
pose into LiBH4 and LiZn2(BH4)5 within this tempera- 
ture range. 
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FIG. 3: (Color online) Formation enegies AF^||, AF^^^, 
AFfooK, and AF^f^ (in units of kJ/mol cation) of Li/Zn 
double-cation borohydrides shown vs. x, the Zn composi- 
tion of the compounds. For LiZn(BH4)3, x = 1/2 while for 
LiZn2(BH4)5, x = 2/3. The structures studied are indicated 
by symbols with the space groups nearby. Diamonds/circles 
are used for theoretically-predicted/experimentally-reported 
phases. Open circle represents the unrefined Cmca structure 
of LiZn2(BH4)5. Constructed from the thermodynamically 
most stable structures, the red-solid lines form a convex hull, 
above which all structures are unstable. 



Figs. 3(a)-(c) demonstrate that at temperatures up 
to lOOK, lattice vibrations have essentially no effect on 
the energetic ordering. Therefore, it is expected that the 
formation energy AF^p^ determined from F^ft ^^^ ^^ 
used to establish more accurately the stability of these 
structures. Fig. 3(d) shows that by considering the vdW 
interactions, the energetic ordering of the /4m2 and the 
Cmca phases of LiZn2(BH4)5 was reverted. In particu- 
lar, the Cmca phase is thermodynamically more stable 
than the IAmi2 phase by c^ 18kJ/mol cation while the 
PI phase of LiZn(BH4)3 is, again, unstable and would de- 
compose into LiBH4 and LiZn2(BH4)5. The instability of 
the PI phase, as shown in Fig. 3(a)-(d), may be the rea- 
son why attempts to experimentally identify LiZn(BH4)3 
have failed. ^^'1^ 

We suggest that the significant role of the vdW interac- 
tions in determining the energetic ordering of the Cmica 
and the /4m2 structures can be traced back to the dif- 
ference in the number of 3D frameworks between them. 




P27c structure 
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FIG. 4: (Color online) Two low-energy structures of 
NaZn(BH4)3: P2i/c (left) and C2/c (right). Sodium and 
zinc and atoms are shown by yellow and gray spheres while 
[BH4]~ complex anions are shown by green tetrahedra. 



Because there are no bonds between the two 3D frame- 
works of the Cmtca structure^^^^ while the /4m2 struc- 
ture has only one 3D framework, they are almost similar 
in their energies at the PBE level. However, in the vdW- 
DF2 calculations, long-ranged interactions between the 
two frameworks of the Cm^ca structure can be captured, 
lowering its energy compared to that of the IAmi2 struc- 
ture. This behavior is similar to that reported by Ref. 
25 where it was argued that the PBE functional artifi- 
cially stabilizes structures of Mg(BH4)2 with unusually 
low densities, and a treatment within non-local density 
functionals such as vdW-DF would be a solution to the 
problem. In this work, we show that to access the stabil- 
ity of LiZn2(BH4)5 at low temperatures, it is obviously 
desirable to properly incorporate the vdW interactions by 
using a non-local density functional such as vdW-DF2. 
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C2/C structure, NaZn(BH4)3 
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FIG. 5: (Color online) Density of phonon states of (a) the 
C2/c structure of NaZn(BH4)3 and (b) the /4m2 structure of 
NaZn2(BH4)5. 
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FIG. 6: (Color online) Two low-energy phases of 
NaZn2(BH4)5: P2i/c (left) and /4m2 (right). Sodium and 
zinc and atoms are shown by yellow and gray spheres while 
[BH4]~ complex anions are shown by green tetrahedra. 



B. Sodium/zinc borohydrides 

For NaZn(BH4)3, the experimentally observed P2i/c 
structure was theoretically predicted^^ to be thermody- 
naniically unstable at OK. A new PI structure was then 
proposed^^ to be lower in energy than the P2i/c struc- 
ture. Using the MHM, we discovered a new monoclinic 
C2/c structure (no. 15), which is lower than the PI 
structure by ^ 3kJ/mol cation. Our further results cal- 
culated with the PEE functional indicate that C2/c is 
the thermodynamically most stable structure at temper- 
atures up to 100 K. The C2/c phase, of which the density 
is p = 1.00 g X cm~^, is slightly lighter than the P2i/c 
phase with p = 1.28 g x cm~^. As an illustration, the 
geometries of the P2i/c and C2/c phases are shown in 
Figure 4. The phonon spectra of these phases were also 
calculated, indicating that both the proposed (PI and 
C2/c) phases are dynamically stable. Fro more informa- 
tion, we show in Fig. 5(a) the density of phonon states of 
the C2/c phase while the density of phonon states of the 
PI phase can be found in the Supplemental Material. ^^ 

NaZn2(BH4)5 was experimentally determined^^ to 
crystallize in a P2i/c phase. The lowest-energy struc- 
ture we discovered by the MHM runs for NaZn2(BH4)5 
belongs to the /4m2 space group, which is thermody- 
namically less stable by just 0.2 kJ/mol cation than 
the P2i/c phase. The IAm2 structure for NaZn2(BH4)5 
can also be obtained by relaxing the /4m2 structure of 
LiZn2(BH4)5 after replacing the Li cations by Na cations. 
Similar to the case of LiZn2(BH4)5, the I4:m2 consists 
of one 3D framework while the P2i/c structure is com- 
posed of two inter-penetrated frameworks. ^^ The I4m2 
structure is therefore a low-density phase with the den- 
sity p = 0.54 g X cm~^, roughly one half of the den- 
sity p = 1.15 g X cm~^ of the P2i/c structure. A geo- 
metrical illustration for the P2i/c and IAm2 phases of 
NaZn2(BH4)2 is given in Fig. 6. The density of phonon 
states of the IAm2 phase of NaZn2(BH4)5, as shown in 
Figure 5(b), indicates that this phase is dynamically sta- 
ble. 

Based on the formation energies, the thermodynamic 
stability of the structures reported for both NaZn(BH4)3 
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FIG. 7: (Color online) Formation energies AP^^|, APq^k ^, 
APfooK, and AP^f.^ (in units of kJ/mol cation) of Na/Zn 
double-cation borohydrides at OK and lOOK shown vs. Zn 
composition x. The structures studied are indicated by 
symbols with the space groups nearby. Diamonds/circles 
are used for theoretically-predicted/experimentally-reported 
phases. Constructed from the thermodynamically most sta- 
ble structures, the red-solid lines form a convex hull, above 
which all structures are unstable. 



and NaZn2(BH4)5 was examined and presented in Fig. 
7. With the PBE functional, the P2i/c and /4m2 
phases of NaZn2(BH4)5 are almost degenerate, being 
the two lowest-energy structures of this compound. For 
NaZn(BH4)3, C2/c is the lowest-energy structure. At 
variance from Ref. 18, the formation energy of the P2i/c 
phase of NaZn(BH4)3 is always negative, indicating that 
this phase is stable and would not decompose into NaBH4 
and Zn(BH4)2. In addition, both the P2i/c and PI 
structures were found to be unstable with respect to the 
decomposition into NaBH4 and NaZn2(BH4)5. The C2/c 
phase, on the other hand, was found to be stable. 

Similar to the case of Li/Zn borohydrides, the vdW 
interactions captured by vdW-DF2 calculations play an 
important role as it lifts the "degeneracy" of the P2i/c 
and the Mm2 structures of NaZn2(BH4)5. In particular, 
the formation energy of the IAm2 structure is slightly 
positive while P2i/c is the thermodynamically most sta- 
ble structure with a strongly negative formation energy. 
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FIG. 8: (Color online) Two low-temperature phases of 
KZn(BH4)3, the R3 phase (left) and PGs/m phase (right). 
Potassium and zinc atoms are shown by purple and gray 
spheres while green tetrahedra represent the [BH4]~ anions. 



To explain this behavior, we note that the IAm2 structure 
consists of one 3D framework while the P2i/c structure 
is composed by two inter-penetrated frameworks that are 
not linked by bonds. The vdW interaction also reverts 
the energetic ordering of the P2i/c and the PI struc- 
tures of NaZn(BH4)3. Moreover, the C2/c phase is again 
the lowest-energy structure of NaZn(BH4)3 and is sta- 
ble with respect to the decomposition into NaBH4 and 
NaZn2(BH4)5. This result shows that from the theoret- 
ical point of view, the existence of the C2/c phase of 
NaZn(BH4)3 is possible. 



C. Potassium/zinc borohydrides 

At 100 K, the experimentally synthesized crystalline 
KZn(BH4)3 sample was determined^^ to be in a rhombo- 
hedral R?> phase. A related compound, KZn5(BH4)5, on 
the other hand, has not been experimentally identified. 
Therefore, it would be interesting to explore, in a similar 
manner to the analysis performed for LiZn(BH4)3, if this 
hypothetical compound is stable or not. Our MHM runs 
for KZn(BH4)3 predicted a hexagonal P6s/m structure 
(no. 176) to be lowest in the DFT energy with the PBE 
functional ^^ft • Compared to the R3 structure, the 
P6s/m structure is thermodynamically more stable by 
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FIG. 9: (Color online) The predicted P2 structure of 
KZn2(BH4)5- Potassium and zinc atoms are shown by purple 
and gray spheres while green tetrahedra represent the [BH4] ~ 
anions. 
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FIG. 10: (Color online) Density of phonon states of (a) the 
PGs/m structure of KZn(BH4)3 and (b) the P2 structure of 
KZn2(BH4)5. 



c^ 3.5 kJ/mol cation. The density of the P6s/m structure 
is p = 1.05 g X cm~^, slightly smaller than that of the R3 
structure with p = 1.15 g x cm~^. For KZn2(BH4)5, a 
monoclinic P2 structure (no. 3) with p = 1.09 g x cm~^ 
was predicted to be quite low in the DFT energy ^^^ft 
(see Table II). The geometries of the P6s/m and R3 
structures for KZn(BH4)3 are shown in Fig. 8 while the 
geometry of the P2 structure of KZn2(BH4)5 is shown in 
Fig. 9. Both the P6s/m phase of KZn(BH4)3 and the P2 
phase of KZn2(BH4)5 were predicted to be dynamically 
stable, as indicated by Figs. 10(a) and 10(b) for their 
densities of phonon states. 

The thermodynamic stability of the R3 and P6s/m 
structures of KZn(BH4)3 and the P2 structure of 
KZn2(BH4)5 was examined and the results are shown 
in Fig. 11. With the PBE functional, the energies of for- 
mation of these structures were determined to be all neg- 
ative at OK and lOOK, implying that they are stable and 
would not decompose into KBH4 and Zn(BE[4)2 [see Figs. 
ll(a)-(c)]. Both the R3 and the P6s/m structures of 
KZn(BH4)3 are located below the line connecting KBH4 
and the P2 phase of KZn2(BH4)5, thus they are also sta- 
ble with respect to the decomposition into KBH4 and 
KZn2(BH4)5. Similarly, the P2 phase of KZn2(BH4)5 
is stable and would not decompose into KZn(BH4)3 and 
Zn(BH4)2. Fig. 11(d) indicates that the inclusion of 
the vdW interactions reverts again the energetic order- 
ing of the R3 and the PQ'^/m structures of KZn(BH4)3. 
Consequently, the experimentally reported R3 structure 
is found to be the thermodynamically most stable phase 
of KZn(BE[4)3. For KZn2(BH4)5, the formation energy 
of the P2 structure remains strongly negative, indicat- 
ing that this phase is thermodynamically stable. The 
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FIG. 11: (Color online) Formation energies AF^ft, ^F^k^, 
AFfooK, and AF^f.^ (in units of kJ/mol cation) of K/Zn 
double-cation borohydrides at OK and lOOK shown vs. Zn 
composition x. The structures studied are indicated by 
symbols with the space groups nearby. Diamonds/circles 
are used for theoretically-predicted/experimentally-reported 
phases. Constructed from the thermodynamically most sta- 
ble structures, the red-solid lines form a convex hull, above 
which all structures are unstable. 



existence of this hypothetical compound KZn2(BH4)5 is 
therefore theoretically possible. 



IV. CONCLUSIONS 

We performed a systematic study on the thermody- 
namic stability of a series of Li/Zn, Na/Zn, and K/Zn 
double-cation borohydrides at low temperatures. Some 
of them were recently synthesized, others are hypothet- 
ical. Using the minima-hopping method, we discovered 



two new lowest-energy structures for NaZn(BH4)3, which 
was experimentally synthesized and KZn2(BH4)5, a hy- 
pothetical compound. The lowest-energy of these com- 
pounds, which belong to the C2/c and P2 space groups, 
respectively, are stable with respect to the decomposi- 
tion into the related compounds, and therefore their ex- 
istence is predicted. For the hypothetical LiZn(BH4)3 
compound, we reach the same conclusion of a previ- 
ous study^^ that the theoretically proposed PI struc- 
ture is lowest in energy. We also found that this struc- 
ture is unstable and would decompose into LiBH4 and 
LiZn2(BH4)5. Therefore, this instability may be the rea- 
son why attempts to experimentally identify LiZn(BE[4)3 
have failed. On the other hand, the experimentally- 
reported structures of LiZn2(BH4)5, NaZn2(BH4)5, and 
KZn(BH4)3 are also predicted to be the thermodynami- 
cally most stable structures. 

While lattice vibrations were found to play a minor 
role in the energetic ordering of the structures studied, 
van der Waals interactions, on the other hand, seem to 
be critically important. We found that the PBE func- 
tional is not sufficient to capture the non-bonding inter- 
actions between the constituent 3D frameworks of the 
complex structures of LiZn2(BH4)5 and NaZn2(BH4)5. 
Consequently, the energetic ordering of the examined 
low-energy structures determined at the PBE level is 
changed. To overcome this problem, we suggest that 
vdW-DF2, a recently developed non-local density func- 
tional, is suitable to access the stability of the low-energy 
structures of these double-cation borohydrides. 
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y^^vTABLE I: Crystallographic information of the lowest-energy structures 
i^iscovered for LiZn2(BH4)5, NaZn(BH4)3, NaZn2(BH4)5, KZn(BH4)3, 
and KZn2(BH4)5. For each structure, cell parameters, the Wyckoff po- 
tiition (WP) and the coordinates (x, y^ and z) of each atom, are given. 
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I4m2 phase, LiZn2(BH4)5 

a = 6 = 8.396 A, c= 16.926 A 

a = /3 = 7 = 90° 



Atom WP X 



y 



Li 2d 


0.00000 0.50000 


0.75000 


Zn 4e 


0.00000 0.00000 


0.13400 


B 8i 


-0.23449 0.00000 


0.18952 


B 2a 


0.00000 0.00000 


0.00000 


H 16j 


0.37333 0.33797 


0.32746 


H Si 


0.35451 0.00000 


0.14971 


H 8i 


0.12618 0.00000 


-0.03731 


H 8i 


-0.24850 0.00000 


0.26071 



C2/c phase, NaZn(BH4)3 
a = 7.793 A, 6 = 15.723 A, c : 
a = 7 = 90°, /3= 110.418° 



7.680 A 



Atom WP X 



y 



Na 4e 


Zn 4e 


B 8f 


B 4e 


H 8f 


H 8f 


H 8f 


H 8f 


H 8f 


H 8f 



0.00000 -0.05711 0.25000 
0.00000 -0.37109 0.25000 
-0.24773 0.04904 0.33521 
0.00000 -0.23185 0.25000 
0.23169 0.47318 0.30728 
-0.09158 0.07548 0.40084 
-0.31184 0.08602 0.18126 
-0.32512 0.06702 0.44784 
-0.36324 0.30786 0.26099 
-0.46155 0.22504 0.39673 



/4m2 phase, NaZn2(BH4)5 
a = 6 = 8.863 A, c = 17.823 A 
a = /3 = 7 = 90° 



Atom WP X 



y 



Na 2d 


0.00000 0.50000 


0.75000 


Zn 4e 


0.00000 0.00000 


0.12723 


B 8i 


-0.22178 0.00000 


0.18054 


B 2a 


0.00000 0.00000 


0.00000 


H 16j 


-0.38010 -0.34645 


0.33640 


H 8i 


0.33654 0.00000 


0.14357 


H 8i 


0.11954 0.00000 


-0.03543 


H 8i 


-0.23158 0.00000 


0.24830 



P63/m phase, KZn(BH4)3 
a = 6 = 8.361 A, c = 7.820 A 
a = ^ = 90°, 7= 120° 



Atom WP X 



y 



TABLE ] 


[ — continued from previous page 


K 


2b 


0.00000 0.00000 


0.00000 


Zn 


2c 


0.33333 0.66667 


0.25000 


B 


6h 


0.18024 0.36048 


0.25000 


H 


12i 


-0.22693 -0.45387 


-0.11496 


H 


6h 


-0.27308 -0.01277 0.25000 


H 


6h 


-0.27307 -0.26031 


0.25000 


P2 


phase. 


KZn2(BH4)5 




a = 


: 7.146 


A, 6 = 5.130 A, c = 


-- 10.151 A 


a = 


:7 = 90°, /3 = 88.741° 




Atom WP 


' X y 


z 


K 


Ic 


0.50000 -0.28780 0.00000 


Zn 


2e 


-0.04620 0.09929 


-0.28148 


B 


2e 


0.29681 0.25477 


0.17862 


B 


lb 


0.00000 0.16741 


0.50000 


B 


2e 


0.16888 -0.14806 


-0.19612 


H 


2e 


0.31447 0.22295 


0.06063 


H 


2e 


0.12578 0.02371 


0.45644 


H 


2e 


0.14255 -0.21993 


-0.08374 


H 


2e 


0.41437 0.38313 


0.23277 


H 


2e 


-0.07793 0.31402 


0.42239 


H 


2e 


-0.14432 0.37736 


-0.19422 


H 


2e 


0.17621 0.09776 


-0.19369 


H 


2e 


0.29818 0.03408 


0.23206 


H 


2e 


0.30971 -0.21894 


-0.25384 


H 


2e 


0.03511 -0.23848 


-0.26096 



to be continued 
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FIG. 1: (Color online) Simulated XRD patterns of the Cmca 
and /4m2 phases of LiZn2(BH4)5. 



FIG. 4: (Color online) Simulated XRD patterns of the R3 
and PGs/m phases of KZn(BH4)3. 
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FIG. 2: (Color online) Simulated XRD patterns of the P2i/c, 
C2/c, and PI phases of NaZn(BH4)3. 
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FIG. 3: (Color online) Simulated XRD patterns of the P2i/c 
and IAm2 phases of NaZn2(BH4)5. 



FIG. 5: (Color online) Density of phonon states of the PI 
phase proposed for NaZn(BH4)3. 
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